Accuracy of bound-state form factors extracted from dispersive sum rules 
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We discuss the extraction of form factors from three-point sum rules making use of harmonic- 
oscillator model, where we derive the exact expression for the relevant correlator. We determine the 
form factor of the ground state by the standard procedures adopted in the method of sum rules, 
and compare the obtained results with the known exact values. We show that the uncontrollable 
uncertainty in the extracted value of the form factor is typically much larger than that for the decay 
constant. In the example considered, we find the uncontrolled systematic error in the extracted 
form factor to exceed the 10% level. 

PACS numbers: ll.55.Hx, 12.38.Lg, 03.65. Ge 

Keywords: Nonperturbative QCD, hadron properties, QCD sum rules 

1. INTRODUCTION 

A QCD sum-rule calculation of hadron parameters [H, involves two steps: (i) one calculates the operator product 
expansion (OPE) series for a relevant correlator and formulates the sum rule which relates this OPE to the sum over 
hadronic states, and (ii) one extracts the parameters of the ground state by some numerical procedure. Each of these 
steps leads to uncertainties in the final result. 

The first step lies fully within QCD and allows for a rigorous treatment of the uncertainties: the correlator in QCD 
is not known precisely (because of uncertainties in quark masses, condensates, a s , radiative corrections, etc.) but the 
corresponding errors in the correlator may be controlled systematically (at least in principle). 

The second step lies beyond QCD and is more cumbersome: even if several terms of the OPE for the correlator 
were known precisely, the hadronic parameters may be extracted from a sum rule only with limited accuracy - the 
corresponding error has to be treated as a systematic error of the employed method. 

In this Letter, we continue our study of the systematic errors of hadron parameters obtained from dispersive 
sum rules. In an earlier analysis [|| we addressed the determination of the decay constant of the ground state by 
means of the two-point correlator. Here, we consider the extraction of the ground-state form factor from the three- 
point correlator in a quantum-mechanical harmonic-oscillator (HO) potential model. This simple model has strong 
advantages compared to more complicated cases: (i) it enables one to calculate the exact three-point function, and 
thus to generate the OPE to any order, and (ii) the bound-state parameters (masses, wave functions, form factors) are 
known precisely. Therefore, we may apply the standard sum-rule machinery to extract the form factor of the ground 
state and then compare it with the known exact form factor. In this way, we may probe the accuracy and reliability 
of the method. (For a discussion of many aspects of sum rules in quantum mechanics, we refer to [1,11, B 00 

We present an explicit example of the form factor extracted at a specific value of the momentum transfer, for which 
the exact correlator is described by the OPE with better than 1% accuracy and standard sum-rule techniques yield 
a form factor extremely stable in the Borel window. However, the value obtained from the sum-rule analysis differs 
from the exact value by more than 10%. 

We therefore reinforce our previous statement that the standard procedures adopted in the method of sum rules 
do not allow one to obtain rigorous error estimates for the ground-state characteristics. In the case of form factors 
extracted from three-point correlators, the uncontrolled systematic errors may be considerably larger than those found 
in the case of decay constants extracted from two-point sum rules. 

2. HARMONIC-OSCILLATOR MODEL 

We consider a nonrelativistic model Hamiltonian H with a HO interaction potential V(r), r = \r\: 

H = H a + V(r), H a =p 2 /2m, V(r) = mwV/2. (2.1) 
The full Green function G(E) = (H — E)~ x and the free Green function Gq{E) = (H — E)~ x are related by 

G- 1 (E)-Go 1 (E)=V. (2.2) 



The solution G(E) of this relation may be easily found by constructing its expansion in powers of the interaction V: 

G(E) = G (E) - Go(E)VG (E) + ■■■ . (2.3) 

In our HO model, all characteristics of the bound states are easily calculable. For instance, for the ground state (n = 0) 
one finds, with q = \q\, 
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E = -u>, i? = |* (r = 0)| 2 = , F (q) = exp(-q 2 /Amw), (2.4) 
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where the elastic form factor of the ground state, Fo(q), is defined according to 



F (q) = <¥ |J($)|* ) = j d?k^(k)^{k-q) = J d 3 r\^(f)\ 2 e^ r (2.5) 

and the current operator J(q) is given by the kernel 

(f'\J(q)\r) = exp(iqr)6W (r - r"). (2.6) 

3. POLARIZATION OPERATOR 

The polarization operator 

II(T) = {f f = 0| exp(- J ffT)|f i = 0) (3.7) 

is used in the sum-rule approach for the extraction of the wave function at the origin (i.e., of the decay constant) of the 
ground state [l[ . A detailed analysis of the corresponding procedure for the HO model was presented in Q . For the HO 
potential, the analytic expression for n(T) is known Q: 
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The OPE series is the expansion of the exact quantity at small Euclidean time T (or, equivalently, in powers of to): 

4. VERTEX FUNCTION 

The basic quantity for the extraction of the form factor in the method of dispersive sum rules is the correlator of 
three currents • The analogue of this quantity in quantum mechanics is 

T{E 2 ,E 1 ,q) = (r f =d\(H-E 2 )-ij(q)(H-E 1 )- 1 \? i =6), q = \q\, (4.10) 

[with the operator J(q) defined in l[2.6p ] and its double Borel (Laplace) transform under E\ — > t\ and E 2 —> t 2 

r(r 2 ,n,g) = (f f = Q\exp(-HT 2 )J{q)exp(-Hn)\fi = 0). (4.11) 
For large T\ and t 2 the correlator is dominated by the ground state: 

r(r 2 , n ,q) - ^ n=0 (r = 0)\ 2 e - E <>^ + ^ F (q 2 ) + ■■-. (4.12) 
Let us notice the Ward identity which relates the vertex function at zero momentum to the polarization operator: 

r(ri,7a,g = 0) =n(n +r 2 ). (4.13) 
This expression follows directly from the relation 

J(q = 0) = l. (4.14) 




In HO model, we find the exact analytic expression for T(j2,Ti,q) by using the results for the Green function in 
configuration space (f' = 0\G(T)\r) from [J]. For our further investigation, we consider the vertex function for equal 
times t\ = T2 = ^T, which has the following explicit form: 
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T(T,q) = [ ex p tanh — . (4.15) 
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The correlator is a function of two dimcnsionless variables ujT and q 2 /mu>. 

Let us now construct for T(T, q) the analogue of the OPE as used in the method of three-point sum rules in QCD. 
The corresponding procedure consists of two steps: First, we expand T in powers of ui 2 and obtain 

r(T,g) = ^r 2n ( g ,rv 2 ". (4.16) 

Each term in this expansion can be computed from the diagrams depicted in Fig. [TJ This is, however, not the full 
story: In three-point sum rules one works with local condensates and, therefore, has for each T2n, n > 1, a power-series 
expansion in T. To keep the same track, we expand T2 n , n > 1, in powers of T. 

In applications, we keep the terms Fq, . . . ,Tg and omit higher-order terms; we then obtain the power corrections 
by expanding T 2 , T 4 , and T 6 in powers of T retaining terms up to order T a . 

As the result of this procedure, the analogue of the OPE for T takes the form 

roPE(T, q) = T (T,q) + T powcl (T,q), 

q 2 T 
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(4.17) 



We display here only terms up to 0(T 4 ) in r powe r but in calculations retain terms up to 0(T 8 ) and O(co 6 ). These 
terms, as well as higher-order terms, may be easily generated from the exact expression (|4.15[) . 

It should be emphasized that the coefficients of each power of T n in the square brackets of (|4.17|) are polynomials 
in q 2 of order (n — 2). Therefore, if the momentum q increases, one needs to include more and more power corrections 
in order to have a certain accuracy of the truncated OPE series for T. In QCD, this implies the necessity to know and 
include condensates of higher dimensions and restricts the applicability of three-point sum rules to the region of not too 
large q 2 . 

Figure [3 demonstrates the behaviour of the exact correlator and the truncated OPE series as described above for a 
fixed momentum transfer q = 1.5 u. Figures [5h.,b make obvious how the ground-state form factor may be extracted 
from the correlator T(T,q) known numerically (e.g., from the lattice): The correlator is dominated by the ground 
state at large values of T; so one may calculate the T- and g-dependent energy 

E(T,q) = -d T \ogT(T,q), 8 T = ^, (4.18) 

which exhibits a plateau at large T: E(T, q) — » E for any q. Making sure that one has already reached the plateau 
and that the correlator is saturated by the ground state, one obtains the form factor from the relation 

F (q)= lim ±e E ° T T(T,q). (4.19) 
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Fig. 2: (a) The energy E(T,qo) = — 8t log F(T, qo) and (b) the correlator T(T, qo) for go = 1-5 w. Blue (full) line: calculation for 
the exact F(T, qo); red (dashed) line: calculation for the truncated OPE series Fope(T, qo) with corrections up to order (ujT) 8 . 
(c) The accuracy of the truncated OPE for T(T, qo) for go = 1.5 u. The vertical line at go = 1.2 a; is the upper boundary of 
the Borel window, (d) Relative contribution of the ground state to the correlator. The vertical line at go = 0.8 a; is the lower 
boundary of the Borel window. 



Fig. [2J; shows that the truncated OPE of (|4. 17|) provides a good (say, better than 1% accuracy) description of T(T, q) 
in the region loT < 1.2. The contribution of the excited states is still rather large in this region of 10T (see Fig. 
and thus a direct determination of the form factor from a truncated OPE is not possible. The procedures of the 
sum-rule method are aimed at modeling the contribution of higher states to the correlator and at obtaining in this 
way the ground-state form factor. 



5. SUM RULE 



The sum rule is merely an expression of equality of the correlator calculated in the "quark" and in the hadron basis: 

R e- EoT F {q) + r cxcitcd (T, q) = T + r power (T, q). (5.20) 

The quantity Tq describes the free propagation and does not depend on the interaction. It may be written as a double 
spectral representation [9(: 

r (T, q) = J dz 1 dz 2 e-? ZlT e-? Z2T A (z 1 ,z 2 ,q), A {z 1 ,z 2 , q) = -^-8 (( Zl + z 2 - q 2 ) 2 - i Zl z 2 < 0) . (5.21) 

Making use of the standard assumption that the contribution of the ground state is dual to the (rectangular) region of 
small values of Z\ and z 2 , we obtain the relation: 

z et l(T,q) z ai{ (T,q) 

R e- EoT F (q)= J dzi J dz 2 e-i ZlT e-^ Z2T A {z 1 ,z 2 ,q)+r powcr (T,q). (5.22) 
o o 

The above expression is exact if we use the exact T- and q-dependent effective continuum threshold, which cannot be 
calculated from the knowledge of only the OPE but can, of course, be reconstructed in our HO model, since we know 
the exact form factor (for details, see Q). Moreover, Eq. (|5.22[) may be even understood as the definition of the exact 
effective continuum threshold, if one makes use of the exact hadron parameters on the l.h.s. Therefore, this sum rule 



5 



alone is not predictive. The form factor (as well as any other parameter) of the ground state may be obtained in the 
method of sum rules only if one imposes an independent criterion to fix the effective continuum threshold. It should 
be, however, understood that this procedure is essentially hand-made and does not arise from the underlying theory. 

The standard assumption is to approximate z e s(T,q) by a T-independent quantity, i.e., to replace it according to 
z e s (T, q) — ► z c . The quantity z c either may be chosen as a q-independent constant or may be adjusted for any value of 
q separately. 

We will now provide an example for the extraction of the form factor from sum rules where all the standard criteria 
point to a very accurate determination of the form factor; the actual error, however, turns out to be much larger. 

Let us consider the extraction of the form factor at qo = 1.5 a;. This specific value is chosen on purpose: for this 
momentum transfer the sum rule for the form factor turns out to be most stable in the Borel window. 
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Fig. 3: (a) The energy of the cut correlator E(T, qo,z c ) = —8t log T(T, qo, z c ) and (b) the sum-rule estimate for the form factor 
at qo = 1.5 u for z c = 2.42a;. (c) The energy of the two-point cut correlator E(T,z c ) = — dr log I1(T, z c ) and (d) the sum-rule 
estimate for the parameter Ro for z c = 2.454oj. The shaded rectangular areas indicate the regions of stability with 1% accuracy. 



First, let us determine the "fiducial" interval (or "window") of T by the following two requirements: (i) The 
truncated OPE gives an approximation to the exact T with, say, an accuracy better than 1%. This yields loT < 1.2. 
(ii) The ground state gives a sizeable contribution of, say, more than 50% to the correlator. This leads to 0.8 < u>T. 
So the "window" where we will work to extract the ground-state form factor is 0.8 < uoT < 1.2. 

Next, we need to impose a criterion in order to fix z c . A widely used procedure is the following (8J: One calculates 

E(T,q ,z c ) ee -9 T lo g r(T,g ,z c ), (5.23) 

which depends on T because of the approximation z e s(T,qo) — > z c (for details, consult 0]). Then, one determines z c 
such that the function E(T, qo, z c ) has a horizontal tangent E = Eq, see Fig. This gives z c = 2.42w, which is used 
to calculate the form factor via Eq. (I5.22|) by the replacement z D g(T, qo) — > z c . Our results are shown in Fig. [3}d. The 
form factor is perfectly flat in the Borel window but nevertheless turns out to be by more than 10% lower than the 
known true value. For comparison, we present in Figs. [3J;,d the corresponding plots for R extracted from the sum 
rule for I1(T) from 0]. Clearly, the general picture in both cases is similar but the deviation from the exact result is 
much greater for the form factor than for the decay constant. 

Let us emphasize a rather dangerous point: (i) a perfect description of T(T, qo) with better than 1% accuracy in the 
"window" , (ii) the deviation of E(T, qo, z c ) from Eq at the level of only 1%, and (iii) a very good stability of F(T, q) 
with better than 1% in the full "window" lead to an error of more than 10% in the extracted value of F(qo)\ Clearly, 
this error could not be guessed on the basis of the other numbers obtained: the full picture mimics a very accurate 
extraction of the form factor, which is, however, certainly not the case. 
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6. CONCLUSIONS 



Let us summarize the lessons one should learn from our analysis: 

1. The knowledge of the correlator in a limited range of relatively small Euclidean times (that is, large Borel masses) 
is not sufficient for the determination of the ground-state parameters. As a consequence, a sum-rule extraction of 
the ground-state parameters without knowing the contribution of the hadronic continuum suffers from uncontrolled 
systematic uncertainties. 

2. Modeling the hadron continuum by a Borel-parameter-independent effective continuum threshold z c allows one to 
fix this quantity z c by, e.g., requiring the average energy E(T) to be close to E in the Borel window. In this case, 
however, the error of the extracted ground-state parameter turns out to be typically much larger than 

(i) the error of the description of the exact correlator by the truncated OPE and 

(ii) the variation of the bound-state parameter in the Borel window. 

3. It is important to realize that the Borel stability of the extracted ground-state parameter — the standard criterion 
that is believed to control both the reliability and the accuracy of the extracted ground-state parameter — does not 
in fact guarantee the extraction of its true value. 

4. The adopted standard procedures for estimating the errors of the extracted bound-state parameters do not allow 
one to provide realistic error estimates. 

The impossibility to control, at present, the systematic errors of the extracted hadron parameters is the weak 
feature of the sum-rule method and an obstacle for using the results from QCD sum rules for precision physics, such 
as electroweak physics. 

Finally, we would like to comment on the obtained quantitative estimates. In HO model, the ground state is well 
separated from the first excitation, which contributes to the correlator, by a large gap of 2ui. This makes the HO 
model a very favourable case for the application of sum rules. Whether or not a comparable accuracy may be achieved 
in QCD, where this feature is absent, is questionable. 
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